Radio Frequency Spectroscopy of Trapped Fermi Gases with Population Imbalance 
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Motivated by recent experiments, we address, in a fully self consistent fashion, the behavior and evolution of 
radio frequency (RF) spectra as temperature and polarization are varied in population imbalanced Fermi gases. 
We discuss a series of scenarios for the experimentally observed zero temperature pseudogap phase and show 
how present and future RF experiments may help in its elucidation. We conclude that the MIT experiments at 
the lowest T may well reflect ground state properties, but take issue with their claim that the pairing gap survives 
up to temperatures of the order of the degeneracy temperature Tf at unitarity. 



PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 

The field of ultra-cold Fermi gases undergoing BCS-BEC 
crossover is particularly exciting because these superfluids 
exhibit a rather novel form of fermionic superfluidity: pair- 
ing begins at temperature T* while condensation take place 
at a significantly lower temperature T c . Concomitantly, the 
normal state fermionic spectrum exhibits an excitation gap or 
"pseudogap" ASS]. Experiments from MIT 10] on popula- 
tion imbalanced Fermi gases simultaneously observe vortices 
and condensate fractions to estimate T c . Combined with RF 
spectroscopy [5] they have thereby established that a pairing 
gap is indeed visible above T c . An even more striking claim 
from this group [5] is that in a highly polarized gas, one finds a 
ground state in which T* ^ while T c = 0. Throughout this 
paper we refer to this state as "a zero temperature pseudogap 
phase". 

It is the goal of the present paper to address these RF spec- 
troscopy experiments on spin imbalanced J3, HI 0, [Ml unitary 
Fermi gases. Our aim is to show theoretically how the self 
consistently calculated RF spectra evolve as temperature and 
polarization are varied thereby accessing the various phases 
which have been contemplated 0, HI fioll . With (roughly) 
decreasing T, these correspond to a Fermi gas, a pseudogap 
phase, a polarized superfluid (or "Sarma" state) and a phase 
separated state. The latter is the ground state for all but possi- 
bly the highest polarizations at unitarity An additional 
aim is to discuss a series of (four) scenarios for the important 
zero temperature pseudogap phase and show how present and 
future RF experiments may help to clarify what is going on. In 
examining one of these we conclude that the lowest T in the 
MIT experiments may, as argued, reflect ground state proper- 
ties. Nevertheless we take issue with another claim from this 
group- that the pairing gap survives up to temperatures of the 
order of the degeneracy temperature Tp. This high value for 
T* at unitarity appears inconsistent with all existing crossover 
theories, and with earlier thermodynamical indications (3J. 

Because it includes the pseudogap in a self-consistent fash- 
ion, we use a theoretical approach J2, [H E2] to BCS-BEC 
crossover which appears to be uniquely positioned to address 
RF experiments lfl3i [l4ll on polarized gases J^, [l5t] at gen- 
eral temperatures T. Presuming, as we do here, the standard 
BCS-Leggett ground state, then Bogoliubov-de Gennes (BdG) 
based schemes are also applicable [16] but only to strictly 
T = 0, where the non-condensed pairs of central interest here, 
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do not enter. Essentially all other schemes in the literature, in- 
spired by the Nozieres-Schmitt-Rink approach 11711 . contain a 
problematic inconsistency in their incorporation of these pseu- 
dogap effects. They presume that [in the fermionic dispersion 
relation E-^ = v/sT+ A 2 (T)] the pairing gap A vanishes at 
and above T c . 

In this paper we use the standard one channel grand canon- 
ical Hamiltonian H — [i\Ni — H2N2 which describes pairing 
between states 1) and 2) and for definiteness take state 1) as 
majority and state |2) as minority, unless indicated otherwise. 
We additionally ignore the interaction between state |3) and 
states |1) and |2), since mean-field energy shifts associated 
with the interaction between |1) and |2) and between |1) and 
3) nearly cancel each other, as observed experimentally. Thus 
state 1 3) is associated with a noninteracting gas. In addition, 
there is a transfer matrix element T k p from |2) to |3) given 

by Ht — Ylk pC^ k .p c 3 p c 2 k + h.c). For plane wave states, 
?k,p = T5(q L + k - p)<5(w kp -wj. Here q L « and u> L 
are the momentum and energy of the RF laser field, and Uk P 
is the energy difference between the initial and final state. It 
should be stressed that unlike conventional quasi-particle tun- 
neling, here one requires not only conservation of energy but 
also conservation of momentum. 

The RF current is defined as I = (N3} = i([H, N3]), Using 
standard linear response theory one finds 



I = 2T 2 lm[X ret {-LU L + ^ 3 - /i 2 )], 
X(iuj n ) = T G 3 (k, %v m ) G 2 (k + q L , iv v 



where ^3 is the chemical potential of |3) and ^23 is the en- 
ergy splitting between |3) and |2). After Matsubara summa- 
tion and using ^(k, v) = 2n8(v — (et + ^23 — M3)) as we U 
as ^2(k, v) = — 2ImC?2(k, v + i0 + ) to rewrite the spectral 
functions for states |3) and |2), respectively, we have 

f 2 

k 

x [/(fik -ui-fj,)- /(e k + ^23 - M3)] , (2) 

where u> = lot, —^23 is defined to be the RF detuning and f(x) 
is the Fermi distribution function. In the above equations the 
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retarded response function X ret (oj) — X(iu> n — » a; + z0 + ), 
and we have expressed the linear response kernel X in terms 
of single particle Green's functions. We define u> n and v m as 
even and odd Matsubara frequencies, respectively and G2 is 
the fully dressed Greens function for the state 2 spins. (We 
use the convention h = ks = Fh 

In our T-matrix formalism lflll[T2ll . G2(k, v) contains two 
self-energy contributions deriving from condensed Cooper 
pairs (S sc ) as well as from finite momentum pairs (S pg ), The 
latter represent pseudogap effects which first appeared in the 
spectral function, in Ref. [18]. We have S = T, pg + E sc , 

A 2 A 2 

where S„„(k, v) = , pg , . and E Sf .(k, v) — ° c . Here 
A sc is the superfluid order parameter, and 7 ^ is associated 
with the life time effects of noncondensed pairs. The resulting 
spectral function can readily be computed as 

A(kv)= 2A 2 g7 (^ + gk) 2 

21 ' ' (y + £k) 2 (* 2 - ^) 2 +7 2 (^ 2 -£-A2 c ) 2 ' 

(3) 

Here £ k ,i = £k - Mi, £k = £k - M> M = (Ml + A*2)/2, 
ft = (/ii — [1-2)1% and v = v — h. In the quasiparticle dis- 
persion, E u , A 2 (T) = A 2 C (T) + A 2 g (T). The precise value 
of 7, and even its T-dependence is not particularly important, 
as long as it is non-zero at finite T. In practice, we choose 
its value based on the experimental atomic peak width. As 
is consistent with the standard ground state constraints, A pg 
vanishes at T = 0, where all pairs are condensed. Above T c , 
we have Eq. © with A sc = 0. Because the energy level dif- 
ference UJ23 (« 80 MHz) is so large compared to other energy 
scales in the problem, the state |3) is initially empty and thus 
7( £ k + ^23 — 1U3) = in Eq. |0. Once the trap is incorpo- 
rated, Eqs. (O and (01 can then be used to compute the local 
current density J(r, u>) and then to obtain the total net current 
I a (u)) = J d 3 r I(r,oj)n a with a = 1,2. Unless stated other- 
wise the energy unit Tp represents the Fermi temperature for 
the noninteracting unpolarized Fermi gas with the same total 
particle number. 

To treat the trap, we assume a spherically symmetrical har- 
monic oscillator potential V(r) = muj 2 r 2 /2. The density, 
excitation gap and chemical potential will vary along the ra- 
dius. These quantities can be self-consistently determined us- 
ing the local density approximation (LDA). The phase dia- 
gram, representing the stable regimes for phase separation, 
the Sarma phase as well as the normal Fermi gas phases as 
a function of temperature and polarization has been mapped 
out |l7lll0[]. Since it is at the heart of the current experments, 
one must also determine H where pairing occurs without su- 
perfluidity. These non-condensed pair effects (which are gen- 
erally ignored in the literature) are also essential for arriving 
at physical values for T c . Important for the present purposes, 
the phase separated state is not associated with pseudogap ef- 
fects, unlike the Sarma state. The same behavior is mirrored 
in the density profiles. The Sarma phase consists of a super- 
fluid core followed by a correlated "mixed normal" or pseu- 
dogap regime, followed by a Fermi gas in the outer regions 
of the trap. The phase separated state, by contrast has an es- 
sentially unpolarized superfluid core separated from a non- 
correlated normal Fermi gas by a sharp interface. The phase 
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Figure 1 : RF spectra for polarized gases in a harmonic trap at uni- 
tarity and polarization 8 = 0.5, for (a) T/T F = 0.4, (b) 0.25, and 
(c) 0.15, respectively. The insets in (b) and (c) are, respectively, the 
pairing peak position and the energy gap A(T) at the trap center as 
a function of T/T^, 1 ' (in units of the majority Fermi energy E^), 
for S = 0.1 (black), 0.5 (red), and 0.8 (green lines), as labeled. The 
corresponding T C /T F = 0.28 (black), 0.25 (red), and 0.19 (green 
lines), respectively, and the estimated T* can also be read off from 
the insets where the gap vanishes. Here we choose 7 = 0.05. 

boundary is determined [ 9] by the balance of pressure. 

To begin, it is useful to present the prototypical behavior 
for the RF spectra. Quite generally we find that in the phase 
separated state (low T) there is a single pairing peak, whereas 
in the pseudogap phase (higher T) there are two peaks. And 
the Sarma phase (intermediate T) may have either one or two, 
depending on T and 5. Finally, at high T, we have only an 
atomic peak, located precisely at the atomic level separation 
UJ23 = 0. For a range of lower T, the atomic peak persists 
deriving from the effectively noninteracting Fermi gas contri- 
bution at the trap edge; the pairing peak arises from the super- 
fluid or pseudogap region in the trap center. 

Figure Q] presents self consistent numerical results for the 
minority RF spectra at unitarity and at moderate polarizations 
6 = (JV-f — Ni)/N — 0.5. The temperature gradually in- 
creases from the lower to upper panels. We consider lower 
temperatures (by a factor of about 2) to arrive at results which 
are comparable to those in [5]. The two insets show the pair- 
ing peak position and trap center gap as a function of T/Tp\ 
In these two insets, we follow the MIT experiments and use 
the majority component Fermi energy Ey as a unit of energy 
for both temperature and gap. The black, red and green curves 
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Figure 2: Low temperature (T = 0.01) RF spectra at different po- 
larizations (as labeled) for a unitary (left) and noninteracting Fermi 
gas. When 5^0, the solid (black) and dashed (red) curves show the 
result when state |2) is the minority and majority, respectively. 



correspond to three polarizations: <5 = 0.1, 0.5 and 0.8 respec- 
tively. One can see that the higher polarization is associated 
with a smaller peak position and energy gap. We see that the 
magnitudes of the pairing gap are rather comparable to their 
experimental counterparts. As in experiment the pairing gap 
increases with decreasing temperature. The energy scale at 
which it smoothly vanishes can be read off in the insets which 
yield T*. There is no sharp feature at T*, so experimentally 
it cannot be precisely defined. Nevertheless, we see that there 
is a clear separation between the peak location curves for the 
three polarizations. By contrast the experimental data for all 
measured polarizations lie on the same (approximately) uni- 
versal curve, with substantially higher T* (by a factor of 2 or 
so). 

We now turn to a first scenario for elucidating the exotic 
non-superfluid phase at high polarizations [5] by considering 
the possibility IU9I1 that this state is a Fermi gas or liquid. The 
loss of superfluidity would be due to a destabilization (arising 
from more benign Hartree-like corrections) in the competing 
normal Fermi gas phase. This scenario is not compatible with 
a zero temperature pseudogap phase (since the presence of an 
excitation gap for fermions means that it is not in a Fermi gas 
or Fermi liquid state). Nevertheless, this scenario would give 
rise to a single, nearly symmetric RF peak at low tempera- 
tures and high polarizations, similar to that observed experi- 
mentally, albeit associated with an atomic rather than pairing 
peak. 

Figure |2] plots the self consistently determined RF spectra 
at very low temperatures T = O.OlTp in the unitary (left) 
and the noninteracting limit (right panels), assuming state |2) 
is the majority (red dashed) and minority (black solid lines), 
respectively. The top two panels correspond to high polariza- 
tions, S = 0.7 and 0.8. The bottom panel presents a com- 



parison with an unpolarized gas. This low temperature phase 
corresponds to superfluidity in all cases in the left column, 
since that is what is found in our self consistent calculations. 
The two high polarizations correspond to phase separation. In 
the noninteracting gas case (right column), the results are very 
simple. We find, as expected, only atomic peaks in the major- 
ity and minority curves. They are located at precisely the same 
position - at the zero of our frequency scale. Comparing the 
two curves in (a) and (b) with (d) and (e) one sees that with 
future majority spectra there is a simple way to rule out this 
particular Fermi gas scenario. At low T the majority curves in 
(a) and (b) (unlike the minority) have atomic peaks as well as 
pairing peaks. The larger atomic peaks of the majority plots 
are associated with the fact that the majority has a much larger 
noninteracting gas tail in its particle density profile. By con- 
trast for the minority curves on the left, all fermions are paired 
at these low T and we see only a single pairing peak. 

When comparing with existing experiments, it should be 
noted that if the single peak in the zero temperature pseudo- 
gap phase were an atomic peak such as in the calculations of 
Ref. 111911 . there would be a shift in its position (relative to that 
computed here) though probably not large enough to match 
the experimental presumed pairing peak. In summary, this 
figure shows that the combined measurement of both majority 
and minority curves can serve to establish whether a single 
peak is coming from paired atoms or noninteracting atoms. In 
this way it can address the scenario which associates the non- 
superfluid state at high polarizations with a Fermi gas phase. 

In Fig. [3] we turn to another possible scenario for the mys- 
terious zero temperature pseudogap phase and calculate the 
self consistent RF spectra within a finite temperature (nor- 
mal) pseudogap phase which arises in the Sarma portion of 
the phase diagram at very high polarizations H. Here we 
consider 5 = 0.95, in order to have polarization and temper- 
atures consistent with the computed phase diagram Js^. The 
temperature gradually increases from T = 0.15Tp to 0.17Tp 
as we go from bottom to top panels. This figure is based on 
the implicit possibility that the purported T = pseudogap 
phase [5] is a finite temperature observation. We thus use RF 
experiments as a type of thermometry and probe whether the 
experimental temperatures are sufficiently low to be in the true 
ground state. Importantly, we see from the figure that a two 
peaked structure is clearly visible at the lowest T of this in- 
termediate temperature scale, 0.152^. It will be even better 
resolved at somewhat lower polarizations, as studied experi- 
mentally. The two peaks start to merge at 0.17Tp, where we 
are left with an atomic peak only. The observation of two 
peaks in this figure, in contrast with experiment, suggests that 
the MIT experiments were conducted at sufficiently low T. 

A third possible scenario for the observed zero temperature 
pseudogap phase follows from BdG-based calculations lfl6ll 
which suggest that the ground state is not phase separated as in 
LDA theories at unitarity, but instead a superfluid with a com- 
plex order parameter - in a Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state 112011 - We have conducted a finite temperature 
study (importantly including noncondensed pairs) of the sim- 
plest such state [21] which suggests that this oscillatory or- 
der parameter phase rapidly becomes unstable with increasing 
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Figure 3: RF spectra at unitarity when |2) is the minority with polar- 
ization S = 0.95, for different temperatures (a) T/T F = 0.17, (b) 
0.16, and (c) 0.15, respectively. Here we take 7 = 0.05. 



temperature. Because it is not sufficiently robust, we argue 



that the FFLO phase is not likely to be a candidate for the ex- 
otic ground state. Indeed, experiments from both groups seem 
to support phase separation flU]. We stress that the phase 
separation that is consistent with the current theory seems to 
be more akin to that in Ref. 16J] than that in Ref. [8] where 
there is very little, if any, pseudogap regime, in either the den- 
sity profiles or the phase diagram. 

As a fourth scenario, we note that the most natural way to 
obtain T c = with T* ^ is associated with phases in which 
there is a frustration of pair mobility which leads to localiza- 
tion of pairs. This state appears in recent theoretical work on 
unitary gases ll22tl in the presence of optical lattices. Never- 
theless, it appears difficult to understand how it can arise from 
Zeeman-like effects, which primarily break pairs apart. Most 
likely (but for very different reasons) this phase has been ob- 
served in high T c superconductors under various perturbations 
where T c , but not T*, is driven to zero 10]. 

In summary, this paper has shown that future RF experi- 
ments are needed to arrive at a more conclusive understand- 
ing of the observed pairing peaks, hopefully, both by reduc- 
ing the unexpectedly high estimates of T* > Tp (in order to 
be consistent with essentially all estimates), and via provid- 
ing majority spectra. The latter can confirm the presence of a 
pairing, as contrasted with, an atomic peak. We cannot rule 
out the possibility that the purported T = pseudogap phase 
has some form of superfluid order. However, if instead a non- 
superfluid but paired ground state is confirmed, it will very 
likely contain some degree of "bosonic" order. 
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